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ABSTRACT 

X-ray absorption spectroscopy provides a potentially powerful tool in determining the metal abun¬ 
dances in various phases of the interstellar medium (ISM). We present a case study of the sight line 
toward 4U 1820-303 (Galactic coordinates l,b = 2°.79, —7°.91 and distance = 7.6 kpc), based on 
Chandra Grating observations. The detection of O i, O ii, O in, O vii, O viii, and Ne ix Ka ab¬ 
sorption lines allows us to measure the atomic column densities of the neutral, warm ionized, and 
hot phases of the ISM through much of the Galactic disk. The hot phase of the ISM accounts for 
about 6% of the total oxygen column density ~ 8 x 10^^ cm“^ along the sight line, with the remainder 
about evenly divided between the neutral and warm ionized phases. By comparing these measure¬ 
ments with the 21 cm hydrogen emission and with the pulsar dispersion measure along the same sight 
line, we estimate the mean o xygen abundances in the ne utral and total ionized phases as 0.3(0.2, 0.6) 
and 2.2(1.1, 3.5) in units of I Anders fc Grevesse I lll989ll solar value (90% confidence intervals). This 
significant oxygen abundance difference is apparently a result of molecule/dust grain destruction and 
recent metal enrichment in the warm ionized and hot phases. We also measure the column density 
of neon from its absorption edge and obtain the Ne/0 ratio of the neutral plus warm ionized gas as 
2.1(1.3, 3.5) solar. Accounting for the expected oxygen contained in molecules and dust grains would 
reduce the Ne/0 ratio by a factor of ~ 1.5. From a joint-analysis of the O vii, O viii, and Ne ix lines, 
we obtain the Ne/0 abundance ratio of the hot phase as 1.4(0.9, 2.1) solar, which is not sensitive to 
the exact temperature distribution assumed in the absorption line modeling. These comparable ISM 
Ne/0 ratios for the hot and cooler gas are thus considerably less than the value (2.85 ± 0.07; Icr) 
recently inferred from corona emission of solar-like stars (Drake & Testa 2005). 

Subject headings: ISM: abundances — X-rays: ISM — X-rays: individual (4U 1820-303) 


1. INTRODUCTION 

The measurement of the ISM metal abundances plays 
a key role in our understanding of the universe. Although 
metals were all produced in stars, the ISM is the deposi¬ 
tory of stellar feedback (e.g., via supernovae). The metal 
abundances of the ISM thus give a unique measure of the 
integrated stellar feedback. The ISM is also a reservoir 
from which stars formed; the similarity or dissimilarity 
of the abundances in the ISM and stars provides im¬ 
portant constraints on the physical processes involved in 
star formation and evolution (e.g., metal diffusion and 
settling). Of course, the abundances are also important 
in the study of various physical and chemical processes 
of the ISM itself (e.g., heating and cooling as well as 
dust grain formation). Furthermore, the interpretation 
of existing observational data often relies heavily on the 
assumption of the abundances in various phases of the 
ISM. Examples are the modeling of hot plasma spec¬ 
tra and the correction for ISM absorption to infer the 
intrinsic spectra of X-ray sources. Therefore, the mea¬ 
surement of the abundances represents a major task of 
astrophysics. 

But a complete accounting of the metal abundances in 
the ISM has been a challenge. The ISM is a heteroge¬ 
neous ensemble of various forms (atomic, molecular, and 
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solid dust grain) and various states (phases): cold ( < 100 
K), warm (^ 8000 K; including warm neutral and warm 
photo-ionized), and hot (^ 10® K; collisionally ionized) 
(e.g., Ferriere 2001 and references therein). For easy of 
reference in the present work, we will not distinguish the 
cold and warm neutral phases and will call their combi¬ 
nation as the “cold” phase. The “warm” phase will, in 
stead, refer only to the above warm ionized phase; the 
total “ionized” phase therefore refers to a combination 
of warm and hot phases. Also these phases are used for 
the atomic gas form only, whereas molecule and solid 
dust grains are denoted as the “compound” form with¬ 
out any phase separation. The abundances of individual 
elements could be significantly different in these forms 
and phases, because of the varying level of depletion and 
molecule/grain destruction as well as recent chemical en¬ 
richment from stellar feedback. The traditional methods 
for abundance measurement in the optical and ultraviolet 
(UV) are sensitive mainly to the cold and warm phases. 
The reliability of such measurements strongly depends 
on how well the physical and ionization conditions are 
modeled (e.g.. Savage & Sembach 1996 and references 
therein). 

Stellar abundance measurements also suffer large un¬ 
certainties, especially for solar-like stars. Recently, so¬ 
lar abundances of light elements such as G, N, O and 
Ne have been revised downward by 25-45 per cent 
llAsplund. Grevesse. fc Sauval Il2005ll from the values of 
Anders & Grevesse (1989; AG89 hereafter). The revised 
values reasonably agree with UV and X-ray measure¬ 
ments of the O abundance in the cold and warm ISM 
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phases (e.g., Sofia & Meyer 2001; Takei et al. 2002; Andre 
et al. 2003; this work). However, this revision has broken 
the abundance accordance with helioseismological mea¬ 
surements and theoretical solar model predictions (e.g., 
Bahcall et al. 2005a). It is argued that the Ne abundance 
in the Sun is poorly determined and that if the Ne/0 ra¬ 
tio is in fact substantially larger (by a factor of > 2.5) 
the models can then be brought back into agreement 
with helioseismological measurements (e.g., Bahcall et 
al. 2005b). Based on the spectroscopy of the O viii, 
Ne IX, and Ne x emission lines in the Chandra HETG- 
ACIS spectra of nearby solar-like stars, IDrake fc Testa I 
1I2OOI estimate that the averaged Ne/0 number ratio is 
0.41 ±0.01 (1 cr), or 2.85 ±0.07 of the AG89 value. How¬ 
ever, two latest independent studie s suggest no such high 
Ne/0 ratio at the solar surface llSchmelz et al. 1 120051: 
lYouugllWil . 

X-ray absorption spectroscopy can, in principle, be the 
ideal tool to measure the metal abundances of the ISM, 
in essentially all forms and phases. The X-ray band¬ 
pass contains almost all the K- and L-transitions of dif¬ 
ferent charge states of the abundant elements from car¬ 
bon to iron. X-ray spectroscopy, less affected by extinc¬ 
tion, also probes larger column densities than is possi¬ 
ble in the optical and UV, and is thus especially use¬ 
ful for measuring the general ISM through much of the 
Galactic disk. With the still limited spectral resolu¬ 
tion and sensitivity of existing X-ray telescopes, how¬ 
ever, only the most abundant elements (like oxygen and 
neon) can produce significant absorption features in the 
spectra of bright background sources such as AGNs and 
X-ray binaries. With the grating instruments aboard 
Chandra and XMM-Newton, several groups have indeed 
studied the abundances of X-ray-absorbing gas, based 
chiefly on measurements of absorption edges, which are 
contributed by the col d and warm phases in both atomic 
and compound form s l|Pa£rels_et_alJ|2003: [Schulz_et_alJ 
2002: iTakei et akl 120021: I.Tuett fc Ghakrabartvl I2003L 


2005:lGunningham. McCrav. fc Snow 112004 lUeda et al. I 


200511 . Most of these studies co uld only get re l ative m etal 
abundances (e.g., Ne/0). But ITakei et al. I II2002I1 and 
ICunningham. McGrav. fc Snow I ll2004ll were also able to 
estimate the absolute abundance (e.g., 0/H) for the sight 
lines to Cyg X-2 and X Persei, by using emission line 
data (e.g, 21 cm) and UV absorption lines. Some success 
has been achieved even in distinguishing the contribu¬ 
tions from the atomic and compound forms through cali¬ 
brating the still uncertai n wavelengths and cross-sections 
of the absorption edges llTakei et al. Il2003l . 

More recently, Yao & Wang (2005; hereafter Paper I) 
and Wang et al. (2005; hereafter Paper H) have demon¬ 
strated the possibility of measuring the abundances of 
the hot ISM phase. By jointly fitting various highly- 
ionized species, significantly detected or not, one can si¬ 
multaneously determine multiple parameters of the mod¬ 
eled absorbers such as the gas temperature, velocity dis¬ 
persion, and element abundance ratio. Even with the 
limited coun t ing st atistics of the LMG X-3 observation, 
iWang et al. I i200,5ll are able to obtain a meaningful lower 
limit to the Ne/0 ratio, which is about the solar value 
of AG89. 

We also note that the detection of O I, O II, and O III 
absorption lines allow for direct measurements of the cold 
and warm oxygen column densities. These lines have 


large oscillator strengths relative to the expected velocity 
dispersions of the ISM phases; there is little fiat part in 
the curve-of-growth (GoG) of such a line. Therefore, one 
can reliably measure the column densities in the typical 
square-root ranges of the GoG. 

Here we capitalize on the high-quality detections of 
multiple absorption lines in the Chandra spectra of 
4U 1820-303 (Futamoto et al. 2004; Paper I) to estimate 
the abundances in the multi-phase ISM of the Galaxy. 
Throughout this paper, we use the AG89 solar abun¬ 
dances as a convenient reference; the number ratios of 
0/H and Ne/0 are 8.5x10“^ and 0.144; in compari¬ 
son, the recently revised solar values are 4.6 x 10“"^ and 
0.151, respectively (Asplund, Grevesse, & Sauval 2005; 
further discussion in § 5). We also assume a collisional 
ionization equilibrium (CIE) in the hot intervening X- 
ray-absorbing gas; this should be quite a good approxi- 
matio n for the oxygen and neon in typical hot ISM condi¬ 
tions llSiitherland fc Uonita, llT99^ . The quoted parame¬ 
ter errors are all at 90% confidence levels unless otherwise 
specified. 

2. 4U 1820-303 AND Chandra DATA CALIBRATION 

4U 1820-303 is an exceptionally good target for the 
X-ray absorption spectroscopy of the ISM. First, this 
low mass X-ray binary (LMXB) is super compact, with 
an orbital period of only 685 s l|Ste11a, et al. Ill 98^ : the 
large ionization parameter in the immediate vicinity of 
the binary greatly reduces the chance for a significant 
local contribution to the observed absorption lines, con¬ 
sistent with their constancy, moderate ionization state, 
and insignificant width/shift (Futamoto et al. 2004; Pa¬ 
per I). Second, the binary resides in the globular cluster 
NGG 6624 {I, b = 2°.79, -7°.91, distance = 7.6±0.4kpc; 
Kuulkers et al. 2003) near the Galactic center (distance 
= 7.62 ± 0.32 kpc; Eisenhauer et al. 2005). Therefore, 
the pathlength of the sight line is well determined and 
samples the entire inner Galactic disk radially to a height 
of ~ 1 kpc off the Galactic plane. Third, the globular 
cluster also contains two radio pulsars, PSR 1820-30A 
and PSR 1820-30B, with essentially the same dispersion 
measure (DM) of 87 cm“^ pc or a total free electron 
colum n density of Ne = 2.7 x 10^° cm“^ llBiggs et al. I 
mni, very useful for an absolute abundance measure¬ 
ment of ionized gas. 

Chandra observed 4U 1820-303 three times: one 
observation used the low-energy transmission grating 
with the high resolution camera (LTEG-HRG; Brinkman 
et al. 2000) at the focal plane for 15 ks, and the 
other two used the high-energy transmission grating 
with the advanced CCD imaging spectrometer (HETG- 
ACIS; Canizares et al. 2005) for 9.7 and 10.9 ks, respec¬ 
tively. The LETG-HRC observation, particularly sen¬ 
sitive to O VII, has been reported by iFiitamoto et al. I 
|[200l . while the HETG-ACIS observations, sensitive to 
Ne IX, have been studied together with other targets by 
Juett, Schulz, & Chakrabarty (2004) and in Paper I. 
These studies have focused on either the global distri¬ 
bution of the intervening hot gas (Futamoto et al. 2004; 

oxygen absorption edge structure 
Here we combine all these obser¬ 
vations of 4U 1820-303, improving both the energy cov¬ 
erage and the counting statistics, to measure oxygen and 
neon abundances in the cold, warm, and hot phases. 


Paper I) or on th e 
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We reprocess the three observations, using CIAO 3.2.1 
and the calibration database CALDB 3.1.0^. We use 
CIAO thread tgextract to extract the source spectra. We 
rebin the MEG spectra to a bin size of 0.0125 A so that 
they can be easily combined with the LETG spectrum. 
For each observation, we calculate the response matri¬ 
ces files (RMFs) and effective area files (ARFs) for the 
positive- and negative-grating arms by running scripts 
mkgrmf and mkgarf. For the HETG-ACIS observations, 
the effective area of the MEG is ~ 5 times larger than 
that of the high energy grating (HEG) in the wavelength 
band where our diagnostic absorption lines are located 
(e.g., longward of ISA; Tabled). The effective area of 
the 1st orders is also at least one order of magnitude 
larger than that of the higher orders (> 2). Therefore, 
we only utilize the 1st order MEG spectra in this study. 
For the LETG-HRG observation, the detector has little 
intrinsic energy resolution, thus one needs to account for 
not only the local line spread but also the global spec¬ 
tral overlapping of different orders. We first calculate 
the RMFs and ARFs from the 1st to 6th grating or¬ 
ders [assuming the higher (> 6) order contribution to the 
spectrum in our interested wavelength range is negligible; 
see below] and then add these RMFs and ARFs pairs to 
form order-combined response files (RSPs; see Nicastro 
et al. 2005 and Paper 11 for more details). To improve the 
counting statistics, we further co-add the positive- and 
negative-grating MEG spectra and combine them with 
the LETG-HRG spectrum to form the final LETG-I-MEG 
spectrum. The corresponding RMFs and ARFs (RSPs) 
are averaged, channel by channel, with a globally mod¬ 
eled spectrum as the weighting function (similar to what 
is described below for the combined spectrum). We con¬ 
duct various consistency checks for key results, based on 
the spectra from individual observations. 

To account for the spectral order overlap, we model 
our combined spectrum over a broad range of 2-24.5 A. 
As in previous studies (Nicastro et al. 2005 and Paper 
11), we fit the continuum with a combination of models: 
a blackbody (BB) and a broken power-law (BKNPL), 
plus ten broad (cr > 1000 km s“^) Gaussians (3 nega¬ 
tive and 7 positive). These broad Gaussian features are 
used to compensate for both the deviation of the data 
from the simple BB-I-BKNPL combination and the cali¬ 
bration uncertainty in various gap/node regions. In ad¬ 
dition, a foreground absorption with solar metal abun¬ 
dances is also applied. This continuum model gives a 
fit with x^/(i.o./=2448/1784. The best-fit hydrogen col¬ 
umn density is Nh = 2.0 x 10^^ cm“^, consistent with 
the values obtained from the optical extinction and H I 
21 cm emission measurements (see § 4). This global fit is 
sufficiently accurate for our order overlap modeling. The 
HRC-LETG higher order contributions are relevant at 
wavelengths longward of 10 A and are ~ 7% and 40% at 
15 and 20 A, respectively. Our subsequent spectroscopy 
is based primarily on local spectral data around indi¬ 
vidual absorption edge/lines considered. The spectral 
analysis uses the X-ray spectral analysis software pack¬ 
age XSPEC (version 11.3.1). 


3. ANALYSIS AND RESULTS 
http://cxc.harvard.edu/ciao/download 


The detected O i, O ii, O iii, O vii, O viii, and Ne ix 
Ka absorption lines (see also Juett et al. 2004; Futamoto 
et al. 2004; Paper I) are clearly visible in our combined 
spectrum. We use six narrow negative Gaussian profiles 
to characterize these lines, which improves the spectral fit 
(x^/d.o.f=2042/1768). Tabled summarizes this Gaus¬ 
sian characterization of the absorption lines, including 
the velocity shifts (cz) relative to the adopted rest-frame 
wavelengths (A) and the equivalent widths (EWs) as well 
as the signifi cance of each line (Ay^). W e note that the 
O II Ka line ll.luett fc Ghakrabartv Il200'^ may be signif¬ 
icantly contam inated by a weak oxygen compound line 
at ^ 23.35 A l|Ta.kei et af. Il200'^ . But the wavelength 
and oscillator strength of the line (or complex) are very 
uncertain, and depend on the exact composition of the 
compound. All the observed lines are unresolved; the 
instrument resolution gives an upper limit ^ 500 km s“^ 
to the velocity dispersion {^/2 times the Gaussian width) 
of the O VII Ka line. We also obtain an upper limit to 
the EW of the O vii K/3 line by fixing its centroid and 
linking its velocity dispersion with that of O vii Ka (Ta¬ 
bled)- This Gaussian characterization of the lines is con¬ 
sistent with previous studies (Futamoto et al. 2004; Juett 
et al. 2004; Paper I), although our present results have 
smaller uncertainties because of the improved counting 
statistics in the LETG-I-MEG spectrum. Since all the 
lines, except for O vii K/3, are saturated, the EWs (pre¬ 
sented here only for ease of comparison with similar anal¬ 
yses) could be strongly biased due to the inaccuracy of 
the Gaussian characterization (Paper I) and should thus 
be interpreted with caution. 

Our analyses and results, presented in the following 
sections, are based on the recently constructed absorp¬ 
tion line model ahsline, as detailed in Paper 1. This model 
adopts the accurate Voigt function as the line profile, ac¬ 
counting for line saturation automatically and allowing 
for a joint-fit to multiple absorption lines. The basic 
model parameters are the line centroid Ei, velocity dis¬ 
persion by, and reference ionic column density Nx (e.g., 
A=0 VII or H). In a fit to multiple lines, other param¬ 
eters can be included, such as the temperature T and 
metal abundance(s) Zx- Such a fit can be much more 
powerful than a simple GoG analysis: e.g., parameters 
such as the velocity dispersion can be measured from 
relative line saturations, and statistical uncertainties in 
multiple parameters can be easily determined. Paper II 
includes a detailed comparison of the absline model with 
the GoG analysis, which is more illustrative but less pow¬ 
erful in the present context. 

3.1. Cold and Warm Phases 

We fit the O i, O ii, and O iii Ka lines, indi¬ 
vidually, with the absline model [Fig. ^ (a)]. We 
assume a velocity dispersion by = 62 km s“^, in¬ 
ferred from the Al iii absorption line in the high- 
resolution far-UV spectrum of HD 163522, a star in 
an adjacent direction of (I, 5)=(350°, -9°) and at 
a dis tance of ~ 9 kpc ijSavage. Sembach. fc Massal 
Eing). The fitted column densities are log[Aioi(cm ^)] 
= 17.6(17.3, 17.9), log[Afoii(cm-2)] = 17.4(16.9, 17.6), 
and log[/Voiii(cm“^)] = 17.0(16.5, 17.5). These values 
are not very sensitive to the exact choice of the assumed 
by value, because these X-ray absorption lines are in the 
saturated square-root range of the CoG (§1); a factor of 
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TABLE 1 

Line parameters and Gaussian fit results 



A 



cz 

EW 

line 

(A) 

fij 

Ax^ 

(km s“^) 

(mA) 

O 1 Ka 

23.508 

0.431 

205 

4(-92, 124) 

83(66, 141) 

O 11 Ka 

23.352 

0.424 

36 

l(-287, 106) 

86(52, 173) 

O III Ka 

23.106 

0.141 

28 

251(-47, 550) 

74(41, 109) 

O VII Kq 

21.602 

0.696 

49 

204(89, 328) 

40(30, 53) 

O VII K/3 

18.627 

0.146 

2 

fixed 

< 10 

O VIII Ko 

18.967 

0.277 

34 

-48(-215, 119) 

19(13, 25) 

Ne IX Kq 

13.448 

0.657 

52 

-35(-213, 68) 

9(7, 11) 


Note. — The rest-frame wavelength (A) values of the O I, 
O II, and O III Ka absorption lines are from averages over multiple 
detections; absolute theoretical values are still quite uncertain (see 
Juett et al. 2004; Paper II). The transition oscillator strength fij 
values of these lines are from |Garcia_e^^IJ 12003). The A and fij 
values of the other lines are as in Paper I. is the reduction 
when the corresponding line is included in the spectral fit. 


two variation in by, for example, would cause a change 
< 20% in these column densities, which is smaller than 
the statistical uncertainties. For O i, in particular, we 
find that the column density changes by < 10% when by 
changes from 62 to 0 km s“^. Systematically, however, 
A^oii might be over-estimated because the potential con¬ 
fusion with compound oxygen absorption (e.g., Takei et 
al. 2002; further discussion in §4). 

These column density estimates, based on absorption 
li nes, may be compar ed with the absorption edge results 
of iJuett et n~l ll2C)04D from the same source. They mea¬ 
sured the column density of oxygen only in the atomic 
form, although the atomic edges are embedded within 
similar structures produced by oxygen in the compound 
form. The compound edge structures are at lower ener¬ 
gies and are less well defined than the atomic ones be¬ 
cause of such uncertain factors as grain and molecule 
composition and confusion with variou s absorption lines 
(I.Tiiett et al. I l200^_|Tajcei_et al_|j2 002D . In addition to 
atomic edges. I.Tuett et al. I ll2004|) also focused on the 
identification and calibration of the Ka lines from O i, 
O II, and O iii. Their analysis was based on HETG ob¬ 
servations of seven X-ray binaries, including 4U 1820- 
303. Toward this source, they first obtained an oxy¬ 
gen column density of 1.31(1.17, 1.51)xl0^® cm“^ from 
the edge measurement and then predicted an EW of 
the O I Ka absorption line. This EW is, however, sig¬ 
nificantly higher than the observed value. Similar off¬ 
sets are also apparent for most of the other sight lines. 
They tentatively attributed this difference to their use of 
the Gaussian line function, instead of the more accurate 
Voigt profile. While the latter profile is used in our analy¬ 
sis, the offset for 4U 1820-303 remains; i.e., our measured 
O I column density 4.0(2.0, 7.9) x 10^^ cm“^ is substan¬ 
tially smaller than the oxygen column density from the 
edge measurement. In fact, an offset is expected because 
the oxygen absorption edge includes contributions not 
only from O i, as is traced by its Ka line, but also from 
the oxygen in the compound form and from O ii and 
O III in the warm phase of the ISM. The sum of our 
measured column densities of O i, O ii, and O iii gives 
N(0 H-0 ii-kO iii)=0.75(0.48,1.22) x 10^® cm-^, which 
is in a better agreement with the edge measurement. 


The remaining small discrepancy is likely caused by a 
slight over-estimate of the atomic oxygen column den¬ 
sity from the edge measurement, presumably because of 
ne glecting the compo und contribution, as acknowledged 
bv l.Iuett et al. ](|20^. Furthermore, their edge model as 
plotted seems to overpredict the absorption at the edge 
(< 23 A; see their Fig. 4). 

We can also measure the column density of neon based 
on its absorption edge. Because neon is a noble ele¬ 
ment and is not easily in a compound form, the absorp¬ 
tion edge structure is much simpler than that of oxy¬ 
gen. But, in the HETG observations, the neon edge 
of the MEG negative grating arm data unfortunately 
lands in a CCD gap, and there are also some compli¬ 
cated instrumental features around the edge wavelength 
of the positive arm data. We thus use only the LETG 
spectrum for the neon edge measurement. We fit the 
spectrum in the wavelength range of 2-24.5 A with the 
XSPEC varabs model, assuming the solar abundances, 
except for the neon, which is modeled separately with 
the XSPEC edge model (Eig. |3). We obtain the neon 
edge wavelength Xe =14.28(14.23, 14.35) A and opti¬ 
cal depth te = 8.6(7.0,10.2) x 10“^. Adopting the 
neon absorption cross section of 3.67 x 10 ~^^ cm“^ from 
fBatucihska-Church fc McCammonI lll992D , we derive the 
neon column density as Vns = 2.3(1.9, 2.7) x lO^"^ cm“^. 

By combining these measurements of the oxygen 
and neon column densities, we further estimate the 
Ne/0 ratio in the cold plus warm gas as (Ne/0)cw = 
VNe/A^(0 H-0 II-LO III) =0.3(0.2, 0.5), or 2.1(1.3, 3.5) 
solar. 

3.2. Hot Phase 

We first consider the X-ray absorption by hot gas of a 
single temperature and a single velocity dispersion (ne¬ 
glecting the small ion mass differences between oxygen 
and neon). Adopting O vii as the reference ion. Ta¬ 
ble 0 illustrates, step by step, the dependence of param¬ 
eter constraints on the inclusion of individual absorption 
lines. The logic in the parameter constraints in each step 
is described in the following: 

1. Fitting the O vii Ka absorption line alone pro¬ 
vides an upper limit to by (determined by the instrument 
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Wavelength A (A) 



Wavelength A (A) 



Wavelength A (A) Wavelength A (A) 

Fig. 1.— The oxygen and neon absorption lines in the LETG+MEG spectrum (with a bin size of 0.0125 A), together with the best-fit 
absline model plus the continuum. The ARE and RMF have been applied to the model. The local jd.o.f. are 156/130 (a), 146/154 (b), 
203/153 (c), and 177/153 (d). The wavelength is marked for the 1st order, although the spectra are also contributed by higher orders. The 
O VII K/1 (c), although not significantly detected, is included in the fit. See text for detail. 



Wavelength A (A) 


Fig. 2. — The Ne edge and the Ne IX Kq absorption line in 
the HRG-LETG spectrum, fitted with XSPEC edge plus absline 
models (x^/d.o./. = 150/144). The bin size is 0.0125 A; the ARF 
and RMF have been applied to the model. The rest is the same as 
in Fig.Q] 


resolution) and a corresponding lower bound to Vqvii 
( whereas its upper limit corresponds to ~ 0 km s“^). 

2. Jointly fitting the O vii Ka and K/3 lines gives a 
lower limit to by and therefore an improved upper limit 
to iVovii- Because both lines trace the same ion, an ex¬ 
tremely small by would make the O vii Ka line too sat¬ 


urated to be consistent with the observed O vii K/3/Ka 
constraint. 

3. The inclusion of the O viii line gives a measurement 
of T (assuming a uniform isothermal CIE plasma origin 
for the lines). 

4. Finally, the Ne ix line, together with the measure¬ 
ment of T in Step 3, constrains the neon column density. 
The comparison of the neon and oxygen column densities 
then leads to the estimate of the Ne/0 ratio. 

The line centroid shifts obtained from the final joint- 
fit [Fig.n (b)-(d)] are nearly the same as those obtained 
from the Gaussian model fits (Table ^1 . The parame¬ 
ters (T, by, and Vovii)) though more tightly constrained 
here, are still consistent with those found in the previ¬ 
ous studies within the quoted error ranges (Futamoto et 
al. 2004; Paper I). But the measurement of the Ne/0 ra¬ 
tio here is completely new. Fig. |21 shows the confidence 
contours of the Ne/0 ratio vs. T. 

One may expect variation in physical and chemical 
properties of the hot gas along the sight line to 4U 1820- 
303. The so-called Local Bubble around the Sun, for 
example, is believed to contain hot gas with T ^ 10® K 
(e.g., Kuntz & Snowden 2000), lower than the average 
value along the sight line (Table |2I). Because the Sun is 
not supposed to be at a “privileged” spot in the Galaxy, 
it is natural to assume that much of the Galactic disk 
is occupied by such bubbles, and that somewhat hot¬ 
ter versions of which may be responsible for the ubiq- 
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TABLE 2 

Step-by-step ahsline fits to hot gas absorption lines 


step 

Included line(s) 

by 

(km s“^) 

log AIovii 
(cm~^) 

log T(K) 

Ne/O 

1 

O VII Ka 

< 446 

17.2(16.3, 18.7) 



2 

O VII Ka, K/3 

298(169, 505) 

16.3(l6.1, 16.5) 



3 

O VII Ka, K/3, O VIII Ka 

325(197, 490) 

16.3(16.1, 16.5) 

6.34(6.29, 6.41) 


4 

O VII Ka, KjS, O VIII Ka, Ne IX Ka 

255(165, 369) 

16.3(16.1, 16.5)“ 

6.34(6.29, 6.41) 

1.4(0.9, 2.1)*' 


Note. — “ The equivalent hydrogen, O VIII, and Ne IX column densities are log[V^f(cm ^)]=20.1(19.9, 20.2), 
log[A^OViii(cni“^)] = 16.4(16.2, 16.6), and log[A^]sieix(cm“^)] = 16.0(15.9, 16.1). ^ The Ne/O ratio is in units of the 
AG89 solar number ratio. 



6.25 6.30 6.35 6.40 6.45 

log[T(K)] 


Fig. 3.— The 60%, 90%, and 99% confidence contours of the 
Ne/O ratio and gas temperature, as obtained in the final joint-fit 
of the hot gas X-ray absorption lines with the absline model (Step 
4 in Table 0- 


niton s 3/4 keV backgronnd emission llSnowden et al. I 
[1993). X-ray absorption line spectroscopy of mnltiple 
sight lines has also shown evidence for the presence of 
a thick Galactic dis k with a scale height of 1-2 kpc 
l|Yao fc Wang ll2?)fi^ . Fnrthermore, the Galactic bulge 
is apparently enhanced in X-ray emission (Snowden et 
al. 1997), although whether it arises mostly from un¬ 
resolved point-like sources or truely diffuse hot gas is 
yet to be determined. A preliminary comparison of 
the 4U 1820-303 sight line with those toward LMC X- 
3 (I, b = 273°.58,-32°.08, D = 50 kpc; Paper II) and 
4U 1957-kll (V1408 Aql; l,b = 5I°.31,-9°.33, U = 5.5 
kpc; Yao et al. in preparation) shows that the averaged 
temperature of the hot gas along the sight line through 
the bulge appears hotter by a factor of ^ 1.5. 

Though without a detailed knowledge about the hot 
gas variation, we can still check its effect on our absorp¬ 
tion line measurements. Because the size of a bubble 
similar to the local one is small (dimension ~ 10^ pc), 
its individual contribution (e.g., IVnv TT <3 x 10^® cm“^; 
iFang. Scmbach. fc Canizares fll2003D ) is negligible, com¬ 
pared to the observed total absorption toward 4U 1820- 
303 (Table O. Therefore, we can characterize the varia¬ 
tion within a bubble and/or from one bubble to another 
along a sight line by using a dispersion around a mean 
parameter (e.g., T). This dispersion may also include 
any Galaxy-wide variation. Here we test how the Ne/O 
ratio may be influenced if the isothermal assumption is 


relaxed by adopting a simple extension to a log-Gaussian 
temperature distribution: 

-( logT- logTo)^ 

dN{T) oc e dlog(r), (1) 

where Tq is equivalent to T in the isothermal case and 
is in units of K, while aiogT is the dispersion of logT. 
After revising our absline model to incorporate this dis¬ 
tribution, we repeat the final joint-fit again. This fit 
gives a 90% upper limit of 0.84 to the new parameter 
(TiogT- As shown in the upper panel of Fig. 0] there 
is a strong correlation between logTo and ciogr; i-e., 
the fit prefers a higher centroid temperature Tq when 
the dispersion is broader. This correlation results from 
the observed strong O viii absorption line, which arises 
only in a narrow and relatively high temperature range 
[logT(K)~ 6.35 — 6.50], compared to O vii, which domi¬ 
nates over a much broader and lower temperature range 
(see Fig. 1 in Paper I). The lower panel of Fig. 01 shows 
that the inferred Ne/O ratio is insensitive to aiogT- A 
similar conclusion is also reached in our test with a power 
law temperature distribution. 

4. OXYGEN ABUNDANCES 

In addition to the preceding X-ray oxygen column den¬ 
sity and Ne/O ratio results, we can also estimate the ab¬ 
solute oxygen abundances (0/H). Such estimates, how¬ 
ever, cannot be made for all individual phases separately, 
due to the lack of required measurements. For example, 
we cannot separately estimate the Ne/O ratios in the 
cold and warm phases, because the Ne edge data do not 
allow for individual measurements of the corresponding 
neon column densities. Table 0 summarizes the results 
for various ISM phases and/or their combinations. To es¬ 
timate an absolute oxygen abundance, we need to know 
the hydrogen column density in the same phase and in 
the same form. Because 4U 1820-303 is located above the 
H I disk of the Galaxy, we can use the hydrogen column 
density from the 21 cm emission line survey in the source 
held, N(H I)= 1.5 x 10^^ cm“^ ijPickev X I;Ockma,n I 
1199011 which may be compared to a total cold hydro¬ 
gen column density N(H I)-|-2V(H2)~ 1.9 x 10^^ cm“^ 
inferred from the optic al extinction of the global clus- 
ter E{B — H) ~ 0.32 ijBohlin. Savasfe. & Drakelll978t 

® We use this easily accessible data base, instead of the 
newer, probably more reliable Leiden-Argentine-Bonn HI survey 
IKalberla et al. ll2Q05fl . The difference (~ 15%) between the two 
survey values is much less than the statistical uncertainties in the 
derived parameters. 
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Fig. 4.— The 68%, 90%, and 99% confidence contours for the 
log-Gaussian temperature distribution. 


TABLE 3 

Summary of Results 




Gas Phase 


Parameter 

cold 

warm 

hot 


column density 

O 17.6(17.3, 17.9) 17.5(17.2, 17.8)“ 16.7(16.5, le.S)*” 

17.9(17.7, 18.1)“ 17.6(17.3, 17.8)'' 

H 21.2“ 20.4''’^ 

Ne _ 17.4(17.3, 17.5)“ _ 16.0(15.9, 16.1) 

Abundance® 

O/H 0.3(0.2, 0.6) 2.0(0.8, 3.6) >0.94 

0.5(0.3, 0.9)“ 2.2(1.1, 3.5)'' 

Ne/O _ 2.1(1.3, 3.5)“ _ 1.4(0.9, 2.1) 

Note. — The column densities are all on a logarithmic scale. The 
values for the cold and warm oxygen phases are obtained from the 
fits to individual lines (O I Ka, O II Ka, and O III Ka), whereas the 
values for the hot phase are obtained from a joint-fit to multiple lines 
(O VII Ka, K/3, O VIII Ka, and Ne IX Ko). “ N(0 II+O III), which 
may be over-estimated because of the contamination from a potential 
compound oxygen line that is not accounted for. N(OVII-I-OVIII). 
“ Cold plus warm gas. '' Warm plus hot (ionized) gas. “ The 21 cm 
emission data are used, f The pulsar dispersion measure is used. ® 
All abundances and ratios are listed in units of the AG89 solar values 
(i.e., O/H = 8.5x10“'' and Ne/O = 0.144). 


IKuiilkers et al. 11200,31) . We estimate oxygen abundance 
in the cold phase as (0/H)„ = N{0 I)/iV(H I) = 
2.6(1.3,5.3) X 10-^ or 0.3(0.2, 0.6) AG89 solar. 

For ionized gas, we can use the free electron column 
density Ne from the pulsar DM (§2). The oxygen column 
densities in the warm and hot phases are N{0 IH-O III) 
and AI(OVII-I-OVIII), neglecting O iv - O vi and O ix, 
which are only important in thermally very unstable “in¬ 
termediate” temperatures (T < 10®'® K) and in proba¬ 
bly rare low-density regions with temperatures > 10® ® 


K. We can then define the oxygen abundances in the 
two phases as (O/H)^, = N{0 II-I-O III)/77(I — i)Ne 
and (0/H)/j = N{0 VIH-O VlU)/^r]Ne, where ^ is 
the hot phase fraction of the electrons and rj = 0.84 
accounts for the contribution from helium. Letting 
(0/H);i = a (0/H)u,, we have ^ = 1/(1 + ar^), where 
VN = N{0 Il-fO III)//V(0 VII-t-0 VIII) ~ 8. We ex¬ 
pect that dust grain destruction and chemical enrich¬ 
ment occurs primarily in the hot phase; i.e., the metal 
abundance in the hot phase should be comparable to, 
or higher than, that in the warm phase. Therefore, 
a > 1; specifically, (O/H)^, = 1.8(0.8,3.1) x I0“® for 
Of = 1 and 1.6(0.7,2.8) x 10“® for a = oo. We thus 
adopt 1.7(0.7, 3.1)xl0“® [or 2.0(0.8, 3.6) solar] to ac¬ 
count for this small a-dependent uncertainty. System¬ 
atically, however, this oxygen abundance could be an 
overestimate, because of the potential contamination to 
the O II line by the uncertain compound oxygen line 
as mentioned in § 3. For instance, if the compound 
oxygen column den sity is ~ 0.5/Voi, as suggested by 
iTakei et al. I ll2002D . then (O/H)^; ~ 9.2 x 10“'* and 
6.7 X 10“*' for a = 1 and oo, respectively. On the other 
hand, (0/H)/j oc a cannot be determined without know¬ 
ing a. Nevertheless, we may estimate the mean oxy¬ 
gen abundance in total ionized gas: [A'(OII-I-OIII) + 
/V(OVII-b OVIlLI/fyA'e = 2.2(1.1,3.5) solar, or ~ 1.1 
solar if the potential contamination from the compound 
oxygen is accounted for. 

5. COMPARISONS WITH OTHER ABUNDANCE 
MEASUREMENTS 

Not all X-ray binaries are suitable background sources 
for ISM measurements. High mass X-ray binaries 
(HMXBs) suffer potentially serious stellar contamination 
to absorption features. The observed X-ray-absorbing 
column density of the HMXB C'yg X-1, for exam¬ 
ple, shows an orbit al dependence l|Ebisawa et al. llT996t 
iSchulz et aDl2002|). Even LMXBs could be problem¬ 
atic. There are previous Ne/O measure ments based on 
absorption edge data of several LMXBs llPaerels et al. I 
120011 iJuett fc Chakrabartv 112003112005D. Higher Ne/O 
ratios (up to a factor of several times the solar value) 
are sometimes reported and are suggested to be local 
to the binary systems, which typically show burst be¬ 
haviors and are proposed to have neon-rich white dwarf 
companions. The discovery of the Ne/O ratio variation 
with time in two of such sources supports this scenario 
j.Iuett fc Chakrabartv1l2005HSidoli et al. II2005I). 

Relatively secure X-ray estimates of the ISM oxygen 
and neon abundances in previous 
only for the sight line to Cyg X-2 
With its Galactic position {l,h = 87°.33,—11°.32 and 
distance = 7.2 kpc), the sight line samples gas primarily 
at the solar radius of the Galaxy, providing an interest¬ 
ing comparison with the inner Galactic disk sight line 
to 4U 1820-303. Similar to 4U 1820-303, Cyg X-2 is lo¬ 
cated at a large off-Galactic plane distance, allowing rea¬ 
sonable estimates of the hydrogen column densities along 
the sight line from 21 cm, Ha, and CO emission lines. 
iTakei et al. I ll2002l) estimated that the oxygen abundance 
is ~ 0.47 ± 0.16 solar in the atomic gas and a factor of ~ 
1.5 higher in the atomic plus compound ISM. Their esti¬ 
mates are based on the absorption edges, which are dom¬ 
inated by the cold and warm phases. Although our oxy- 


studies are available 
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gen measurements are based on the absorption lines, we 
can derive an average oxygen abundance for the cold plus 
warm gas along the sight line to 4U 1820-303 (Table 13), 
N{0 I+O II+O III)/[fV(HI) + (l - Or7W]=0.52(0.33, 
0.85) solar. This value is consistent with the oxygen 
abundance estimate along the sight line to Cyg X~2. 

Our neon abundance measurement of the cold plus 
warm gas along the sight line to 4U 1820- 303 is based 
on th e absorption edge, as it is to Cyg X-2 llTakei et al. I 
l 2 ?iol . The neon abundance 0.75(0.55, 0.95) toward this 
latter source is considerably lower than 1.2(1.0, 1.4) in¬ 
ferred here for the sight line to 4U 1820-303, which in¬ 
dicates a moderately metal enhancement in the inner re¬ 
gion compared to the outer region of the Galaxy. The 
Ne/0 ratios along these two sight lines are compatible 
[2.1(1.3, 3.5) vs. 1.6(0.9, 2.3)]. For 4U 1820-303, we also 
have the hot Ne/0 ratio 1.4(0.9, 2.1), comparable to the 
value in the cold plus warm gas, particularly if the ex¬ 
pected oxygen compound contribution is included. The 
hot Ne/0 ratio, if typical in the Galaxy, can be used in a 
joint analysis of oxygen and neon X-ray absorption lines 
observed in other sight lines. 

Now let us compare our inferred ISM Ne/0 ratios with 
those stellar measurements. At first glance, our Ne/0 ra¬ 
tio 2.1(1.3, 3.5)] obtained for the cold plus warm gas 
is consistent with the high Ne/0 ratio as suggested by 
IDrake fc Testa I ll200,5ll . But, including as well the ex¬ 
pected oxygen compound contribution would then lead 
to a reduced Ne/0 ratio (by a factor of ~ 1.5; Takei 
et al. 2002) that would be significantly lower than the 
stellar ratio. Similarly, our inferred hot gas Ne/0 ratio 
(Table|3 is significantly lower (by ~ 3 cr) than the stellar 
value. It is difficult to reconcile this Ne/0 ratio difference 
and the total oxygen abundance similarity between the 
stellar and ISM measurements. There is no obvious way 
to hide neon in the ISM! We suspect that the Ne/0 ratio 
estimate from stellar emission lines, which is generally 
sensitive to modeling assumptions, may not represent 
the true underlying composition of the s tars. The incon¬ 
sisten cy of the stellar measurementej3y_ |Drake_fcTeste I 
120051) with t he recent solar studies lSchmelz^^L^200^ 
lYoung ll200,5l) illustrates such an uncertainty. 

6. SUMMARY 

We have used the Chandra HETG and LETG observa¬ 
tions of the LMXB 4U 1820-303 as a test case to explore 
the potential for metal abundance measurement in differ¬ 
ent atomic phases of the ISM. We have concentrated on 
the measurements of atomic oxygen and neon, as sum¬ 
marized in Table |3 and on comparison with previous 
X-ray studies. Our main results and conclusions are as 
follows: 

1. We have separately measured the column densities 
of O I, O II, and O iii from their heavily saturated Ka 
lines. These measurements are only weakly dependent on 
the exact velocity dispersion assumed. A comparison of 


Anders, E., Grevesse, N. 1989, Geochim. Gosmochim. Acta, 53, 
197 

Andre, M. K., Oliveira, C. M., Howk, J. G., Ferlet, R., Desert, 
J.-M., Hebrard, G., Lacour, S., Lecavelier des Etangs, A., Vidal- 
Madjar, A., & Moos, H. W. 2003, ApJ, 591, 1000 
Asplund, M., Grevesse, N., & Sauval, J. 2005, ASPC, 336, 25 


the O I column density and the 21 cm hydrogen emission 
in the field gives a cold (neutral) oxygen abundance of 
0.3(0.2, 0.6) solar (AG89). The ratio of our measured 
O II plus O III column density to the pulsar DM along 
the same sight line further gives an estimate of the warm 
oxygen abundance as 2.0(0.8, 3.6) solar. 

2. We have constrained the neon column density from 
its absorption edge, giving an abundance of 1.2(1.0, 1.5) 
solar and a Ne/0 ratio of 2.1(1.3, 3.5) solar for the cold 
plus warm gas. The Ne/0 ratio with the inclusion of 
the compound ISM is likely to be a factor of ^ 1.5 lower 

3. We have also measured the oxygen column density 

and Ne/0 ratio in the hot ISM, based on a joint-analysis 
of the detected O vii Ka, O viii Ka, and Ne ix Ka 
absorption lines, together with the non-detection of the 
O VII K/3 line. Assuming that O vii Ka and O viii Ka 
trace all hot gas (i.e., T ~ K), we estimate that 

the hot phase accounts for about 6% of the total oxygen 
column density along the sight line and obtain a Ne/0 
ratio of 1.4(0.9, 2.1) solar, which is insensitive to the 
exact temperature distribution assumed. 

4. Our abundance estimates for the atomic phases, to¬ 
gether with complementary X-ray spectroscopic studies 
of the total abundances, may have strong implications for 
understanding various chemical enrichment and deple¬ 
tion processes both in the ISM and during star formation. 
There is evidence for an enhanced enrichment in the hot 
gas, especially in comparison with the cold phase. The 
Ne/0 ratios obtained in the ISM are significantly smaller 
than the value indicated in t he recent emission lin e mea¬ 
surement of solar-like stars llDrake fc Testa 1120051) . 

The existing X-ray measurements of the ISM abun¬ 
dances (including the present work) are still preliminary, 
subject to various systematic uncertainties, both theo¬ 
retical and observational, which are yet difficult to fully 
quantify. Nevertheless, the outcome of this study demon¬ 
strates the unique potential for a comprehensive charac¬ 
terization of the metal abundances in various ISM forms 
and phases. 
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